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Abstract 



Notwithstanding the long and distinguished history of studies of vibra- 
tional energy relaxation, exactly how it is that high frequency vibrations man- 
age to relax in a liquid remains somewhat of a mystery. Both experimental 
and theoretical approaches seem to say that there is a natural frequency range 
associated with intermolecular motion in liquids, typically spanning no more 
than a few hundred cm -1 . Landau-Teller-like theories explain rather easily 
how a solvent can absorb any vibrational energy within this "band" , but how 
is it that molecules can rid themselves of superfluous vibrational energies sig- 
nificantly in excess of these values? In this paper we develop a theory for such 
processes based on the idea that the crucial liquid motions are those that most 
rapidly modulate the force on the vibrating coordinate — and that by far the 
most important of these motions are those involving what we have called the 
mutual nearest neighbors of the vibrating solute. Specifically, we suggest that 
whenever there is a single solvent molecule sufficiently close to the solute that 
the solvent and solute are each other's nearest neighbors, then the instanta- 
neous scattering dynamics of the solute-solvent pair alone suffices to explain 
the high frequency relaxation. This highly reduced version of the dynam- 
ics has implications for some of the previous theoretical formulations of this 
problem. Previous instantaneous-normal-mode theories allowed us to under- 
stand the origin of a band of liquid frequencies, and even had some success in 
predicting relaxation within this band, but lacking a sensible picture of the 
effects of liquid anharmonicity on dynamics, were completely unable to treat 
higher frequency relaxation. When instantaneous-normal-mode dynamics is 
used to evaluate the instantaneous pair theory, though, we end up with a mul- 
tiphonon picture of the relaxation which is in excellent agreement with the 
exact high-frequency dynamics — suggesting that the critical anharmonicity 
behind the relaxation is not in the complex, underlying liquid dynamics, but 
in the relatively easy-to-understand nonlinear solute-solvent coupling. There 
are implications, as well, for the Independent-Binary-Collision (IBC) theory 
of vibrational relaxation in liquids. The success of the instantaneous-pair ap- 
proach certainly provides a measure of justification for the IBC model's focus 
on few-body dynamics. However, the pair theory neither needs nor supports 
the basic IBC factoring of relaxation rates into many-body and few-body dy- 
namical components — into collision rates and relaxation rates per collision. 
Rather, our results favor taking an instantaneous perspective: the relaxation 
rate is indeed exercise in few-body dynamics, but a different exercise for each 
instantaneous liquid configuration. The many-body features therefore appear 
only in the guise of a purely equilibrium problem, that of finding the liklihood 
of particularly effective solvent arrangements around the solute. All of these 
results are tested numerically on model diatomic solutes dissolved in atomic 
fluids (including the experimentally and theoretically interesting case of I2 
dissolved in Xe). The instantaneous pair theory leads to results in quantita- 
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tive agreement with those obtained from far more laborious exact molecular 
dynamics simulations. 
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I. INTRODUCTION 



The most crucial step in a solution-phase chemical reaction is almost always a thermal 
activation; improbable as it may be, the solvent must find some avenue for concentrating 
enough of its kinetic energy within the vibrations of the reactants to allow the reaction to 
proceed.0! A similar challenge can be posed later on in the process; the ability of the resulting 
products to return any excess vibrational energy back to the solvent is often critical to the 
outcome of the reaction.^! Beyond these issues however, there is a conceptual importance 
to understanding such vibrational energy transfers that goes to the heart of how one thinks 
about the dynamical role of a solvent. If the entire solvent-solute system is regarded as a 
"super molecule" (as might be appropriate for solid-state vibrational relaxation), then the 
question is how a V-V relaxation occurs — that is, which vibrations of the system as a 
whole (which phonons) can serve as repositories for the solute's vibrational energy.@~@ If, on 
the other hand, we consider the solute to be buffeted by occasional collisions with individual 
solvent molecules, a very natural perspective for the gas phase, what we need to understand 
is a V-T processJj So how is it that we should think about vibrational energy relaxation in 
liquids?^ 

A formal answer to this Question arises quite simply from a Fermi's Golden Rule evalu- 
ation of the relaxation rate.ca Within the apparently quite broad reach of low-order time- 
dependent perturbation theory, the rate at which vibrational energy relaxes from a given 
solute mode, Tf" 1 , is proportional to the Fourier component of the autocorrelation function 
of the force on that mode at the mode's own frequency, ujqM For a diatomic molecule of 
reduced mass fi, for example, the classical vibrational energy relaxation rate is predicted to 
obey the Landau-Teller formula,0 

^ = -vW , (i.i) 

where 

fj(Uo) = / dt 7] (t) COS LU t (1.2) 





V (t) = ^ (6F(t) 6F(0)) . (1.3) 

Here ks is Boltzmann's constant, T is the temperature, and by 5F(t) = F(t) — F, the 
fluctuating force on the mode, we mean the hypothetical classical time evolution of the 
solvent that would occur were the mode itself held motionless. 

Leaving aside such issues as whether a purely classical treatment of the solvent can 
adequately represent the influence on what is frequently a rather quantum mechanical vi- 
bration,!^ it seems clear that this approach can usually be counted on to give a numerically 
accurate rendition of classical vibrational energy relaxation.!!! More than that, though, it 
provides a conceptually suggestive rendition because it fits so nicely into the broader con- 
text of relaxation in liquids. Within a Generalized Langevin representation of the problem, 
the function rj(t) can be interpreted as an approximation to the solvent-induced friction felt 
by a vibrating coordinate — and since the exact friction can sometimes be determined by 
simulation, it has been possible both to check and to understand the formalism's accuracyJiHl 
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Still, the fact that numerical success does not necessarily translate into conceptual insight 
becomes clear once we try to get at the central V-V vs. V-T dichotomy. It has long 
been appreciated that any frequency-domain friction f)(u)) can be regarded as having its 
origin in the dynamics of a dense set of harmonic oscillators, oscillators whose frequencies 
span the frequency range of the friction and which are coupled linearly to the coordinate 
of interest .Ho The fact that liquid frictions interpreted in this fashion typically have a 
spectrum spanning a few hundred cm~ l could thus be taken as evidence that there is indeed 
something resembling a phonon band of this width in liquids. It is then but a short step 
to suggest that any vibrational relaxation of modes with frequencies lying within this band 
ought to be considered as occuring via a resonant V-V transfer. 

This argument is actually rather weak as it stands. These harmonic oscillators could, 
in fact, be little more than mathematical constructs lacking any microscopic, molecular 
foundation. What gives the argument a bit more substance is that recent work has shown 
that it is possible to associate these oscillators with the instantaneous normal modes of the 
solutions, oscillators which really are specific, well-defined molecular motions.il~H Indeed, 
perhaps it is fair to say that the basic success of such instantaneous-normal-mode (INM) 
approaches in predicting most of the essential behavior of the vibrational friction lends a 
certain amount of credibility to what would otherwise be fanciful attempts to take liquid 
phonon pictures literally in theories of vibrational relaxation.il 

There is, however, a major problem and it is this problem that forms the topic of this 
paper. Once the vibration to be relaxed lies outside the putative band of the liquid, any 
theory with nothing but linear forces and harmonic oscillators will predict that energy relax- 
ation is completely forbidden.ill Yet, what we might call the high frequency modes of solutes 
(those with frequencies higher than the Debye frequency in the solid-state parlance) really 
do relax, albeit very slowly&0S It is worth noting that the Landau-Teller formula itself 
has no difficulty with this phenomenon; it simply ascribes the relaxation to the presence of 
small but finite high-frequency Fourier components of the force autocorrelation function.0 
The conceptual difficulty is that we have independent evidence about the intrinsic band- 
width of intermolecular vibrations in a liquid, both from theoretical predictions, such as 
instantaneous-normal-mode analysisII~0 and Fourier representations of simulated velocity 
autocorrelation functions^! and from experimental results, such as those from far infrared, 
Rayleigh and Optical Kerr effect spectra.^S None of these estimates give any evidence that 
the band for simple liquids extends significantly beyond a few hundred cm" 1 . A much more 
likely scenario than one that would require every simple liquid to have the kinds of 500-2000 
cm -1 intermolecular vibrations necessary for simple V-V processes is that high frequency 
vibrational energy relaxation is a strongly anharmonic event — and it is this anharmonicity 
we need to understand. 

Having made this point, though, there is an interesting distinction we can make. The 
relaxation of impurity molecules in crystals can often be modeled as a multiphonon process, 
one in which the essential harmonicity of the vibrations (and hence the very idea of phonons) 
is respected, but in which the nonlinearity of the coupling forces allows for V-V transfer to 
overtones of the crystal's intrinsic frequencies. In systems as inherently anharmonic 

as liquids, however, we could easily imagine that the basic idea of independent harmonic 
excitations could be the first casualtyJ£0c3 Thus, without further work, it is far from clear 
which of these two kinds of anharmonicity — the nonlinearity of the forces being driven 
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by the liquid's dynamics or that of the actual dynamics itself — ought to be more critical 
to high-frequency relaxation in liquids. In particular, we are still left with the question of 
whether the V-V paradigm for resonant energy transfer is going to be any better than a 
V-T based perspective would be. 

Complicating this issue is the fact that there is some evidence on the side of molecular 
translation being the primary sink for vibrational energy. That is, there is some support 
for a collisional model for vibrational relaxation. While it is admittedly far from intuitive 
how discrete, uncorrelated collisions could ever occur in a dense a medium as a liquid, 
the so-called Isolated-Binary Collision (IBC) model0El has been remarkably successful in 
rationalizing the thermodynamic state dependence of vibrational population relaxation rates 
based on just such a picture.il'illl The notion behind this model is that these rates should 
factor into a purely few-body scattering piece specifying the rate of energy transfer per 
collision, and a liquid-state piece quantifying the rate at which the solvent collides with the 
solute. The fact that the concept (and thus the rate) of collisions is ill-defined in liquids is 
certainly vexingJ^S but it is hardly out of the question that some other, more sophisticated 
analysis of the collision concept might find that the V-T motif is still a dominant mechanism 
for relaxation. 

The approach we shall take in this paper to getting at these issues is to begin with 
the standard perturbative (Landau-Teller) expression in terms of the force autocorrelation 
function, Eqs. (|1.1| )- (|1.3| ), but then ask from precisely where we should expect the ma- 
jor high frequency components of the correlation function to arise. For any instantaneous 
configuration of the liquid, this analysis will immediately lead us to look at the short-time 
scattering dynamics of the solute and a single crucial solvent partner. Previous work on 
vibrational relaxation has shown that whenever there is a solvent molecule that is nearer 
to the solute than any other solvent molecule is to either member of the solute-solvent 
pair — that is, whenever there is a mutual-nearest-neighbor solute-solvent pair — then 
the high-frequency solvent effects are dominated by the pair, at least within the liquid's 
band. El We will show how this observation can be exploited more generally, transcending 
the limitations to dynamics within the band. Not surprisingly, our theoretical prediction 
for the relaxation rate will then separate into a straightforward few-body dynamical calcu- 
lation and a residual many-body problem, just as with IBC theory. However, unlike the 
IBC approach, the liquid-state portion of our theory will turn out to rest on a simple, well- 
defined, purely equilibrium question. Moreover, our development will make it clear how we 
can make connections between this V-T theory and the harmonic instantaneous dynamics 
of instantaneous-normal-modes, and thereby connect with the V-V energy transfer ideas 
which seem to work so well within the solvent band. 

The remainder of the paper will be organized as follows: Section [Ij will present our 
basic theory, an expression for the full force autocorrelation function of the system in terms 
of a correlation function for single degree of freedom, one whose equations of motion we 
proceed to derive. Section III then explores the connection to the instantaneous-normal- 
mode perspective, showing how we can view the new theory as a natural generalization of 
INM ideas. Section IV provides some numerical examples of how and how well the various 
approaches work, and we conclude in Sec. V with some general comments on what we are 
now in a position to say about the questions posed in this Introduction. 
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II. INSTANTANEOUS PAIR THEORY FOR THE VIBRATIONAL FRICTION 



A. Solute Solvent Coupling versus Solute Solvent Dynamics 

While the central quantity of interest in vibrational relaxation is the force autocorrelation 
function of Eq. ( |1.3|) , it is convenient to look instead at a mathematically equivalent function, 
the "force-velocity" autocorrelation function,H@ 

G F (t) = (F(t) F(0)) = -(k B T) i}{t) , (2.1) 

which lets us write the frequen£y-domain vibrational friction we need in terms of the cosine 
transform of the new functionilllH 

*»> = j^(^f* G 'M~* s iM^)- (2 ' 2) 

The main reason we bother with this conversion is that this alternative version will suggest 
a natural division within our problem. As one might have expected, there will always be 
contributions from the relative solute-solvent velocities — the intrinsic dynamics of the prob- 
lem — but there will be also be separate features dependent directly on the solute-solvent 
distances — and it is these that will be most pertinent to the solute-solvent coupling.il 
Indeed, one of the key aspects of our subsequent development will be the realization that 
these two kinds of information actually play substantially different roles in determining the 
high-frequency behavior. 

We can begin to see how this division occurs by looking at the problem of evaluating 
a general correlation function of the form of Eq. ( p.l|) when F is any specifically solute- 
centered quantity. That is, suppose we look at a correlation function of the form of Eq. 
Q2.1|) with the function F given by a sum of solute-solvent pair "potentials", 

F = E ™(roj) , (2.3) 

j £ solvent 

with the tagged solute labeled as particle 0, the remaining solvent denoted by j = 1,...,N 
and the interparticle distances given by r 0j - = |r — rj|, so that 

F = £ w'(r 0j ) r 0j ■ (2.4) 

j £ solvent 

This general function w(r), something we have referred to elsewhere as the "spectroscopic 
probe potential" ,0 could very well be the force on some portion of the solute (as it is in 
vibrational relaxation), but for the purposes of this argument it could also be the change in 
the solute-solvent potential energies when the solute is promoted from one electronic state 
to another (the situation in time- dependent fluorescence studies of solvation) .e3 

If we substitute Eq. (|2.4| ) into Eq. (|2.1| ), we find that the velocity form of the correlation 
function consists of two kinds of terms, 

G F (t) = V E™W*)] w'[r 0j (0)] r 0j (t) f 0j (0)\ 
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+ ( E ™'K(t)] w'[r ok (0)} r 0j (t) r ok (0) 



(2.5) 



a binary piece, the first term, which considers only the dynamics caused by the motion 
of a single solvent at a time, and a ternary piece, the second term, which involves the 
simultaneous motion of two different solvents. 

These two parts will not contribute equally, for reasons we shall come to presently, but 
the essential ingredients of the correlation function are already visible. When we evaluate 
any of the terms at time t — 0, 

(w'[r 0j (0)} w'[r ok (0)} r Oj (0) r 0fe (0)> = (w'[r 0j (0)} w'[r Ok (0)}) (r 0j (0) r ok (0)) 

= S jk (w>[r Oj (0)} w'[r Oj (0)}) (f 0j (0) r ,(0)) , (2.6) 

we notice that the contribution factors rigorously into a static part involving the coupling, 
w'(r), and a dynamical part involving the velocity, f. This separation suggests, though, 
that this same conceptual distinction between coupling and dynamics ought to continue 
to operate at short times - - and therefore at high frequencies - - despite the fact that 
the expression no longer factors. But, once we make such a distinction, we can simplify 
our expressions dramatically. For one thing, at the highest frequencies we should expect 
coupling to be far more influential than dynamics in determining the correlation function. 
With typical intermolecular potentials and forces, the sharply varying repulsive parts allow 
for tremendous variation in the coupling (and therefore enormous values of w'(r)) just from 
small changes in internuclear distances. The velocities, by contrast, never become all that 
large in magnitude, guided as they are by the equilibrium Boltzmann distribution. 

We can make this same point a bit more quantitatively. If we assume that the t = 
factorization holds, at least roughly, then each term in Gp{t) can be written 
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T jk (t) = (w'[r oj (t)\ w'[r ok (0)} r 0j (t) r ofc (0)) 



« (w'M*)] u/[r„k(0)]) (r 0i (t) r ofc (0)) 
which becomes, on cosine transforming the resulting product, 



T ik ~ V io dUJ ' V ^ \- C<kUJ ~ ^ + C<kU + ^ 

where V(u) and C(u) are the cosine transforms of the single-solvent velocity and coupling 
correlation function, respectively. What then determines the high frequency behavior? The 
intrinsic frequencies of the liquid are largely those of V(u')f^ which cuts off at frequencies 
uo' well below the uo we are interested in. But if u' -C oj, then C{u — to') ~ C(u), so 



2 

Tjk(u)) ~ $jk— 

71 







oo 

duo' V{uj') 



C{u) . (2.7) 



Hence, we see again that the high-frequency dependence of our friction spectrum results 
almost entirely from the high-frequency behavior of the coupling. 

This argument can be taken still further. The dynamic range of the coupling is so large 
that we might expect it to be dominated by the very largest contributions. No more than a 
few of the closest solvents should therefore play much of a role. Indeed, we might venture 
that only the single largest contribution, from the single most strongly coupled solvent is 
important ,EiEI 

G F (t) = (w'[r oj (t)} w'[r Oj (0)} f 0j (t) f oi (0)>. 

=maximally coupled ' \ " / 

What we suggest here is that this latter idea is virtually exact; that the high-frequency 
dynamics of a solute-centered correlation function stems almost entirely from the binary 
dynamics of a single solute-solvent pair — and all we need do is compute the nonlinear, 
anharmonic dynamics of this pair for each instantaneous liquid configuration. 

There are actually a number of kinds of evidence supporting our analysis. An immediate 
consequence of our argument is that the binary terms in Eq. ( |2.5| ) ought to be significantly 
larger than the ternary terms at short times. Indeed, exact molecular dynamics simulations 
have established the overwhelming predominance of the binary contributions to the Gp(t) 
vibrational relaxation correlation functions.il Interestingly, one finds the predicted strong 
preferences for the (presumably) less collective binary part even when the relaxation takes 
place in dipolar solvents, liquids whose intermolecular forces are sufficiently long-ranged that 
the friction r](t) itself ends up with major contributions from ternary terms.0il The reasons 
for this difference between binary and ternary are clear from a short-time perspective. At 
equilibrium (t — 0), the force-velocity version of the vibrational friction is rigorously binary, 
because, as we can see from Eq. ( |2.6| ), the velocities of different solvent molecules are 
rigorously independent. The actual equilibrium friction, by contrast, does not even start 
out as binary at t = 0; the equilibrium structure of the liquid builds in correlations that 
have nothing to do with the subsequent dynamical evolution. But even away from t = 0, the 
wide range of values possible for the couplings serves to explain the observation. The sheer 
unlikelihood of having two different solvents so similar in their distance from the solute that 
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they end up making comparable contributions will also prevent much of a ternary presence 
in the answer. 

We can see somewhat more direct evidence of the two-particle character of the essential 
dynamics by looking in the frequency domain with the aid of INM theory. INM theories 
cannot addresss any of the truly high-frequency dynamics, but if we confine ourselves to 
the highest frequencies inside the INM band, we invariably find that the vibrational friction 
spectrum (p v ib(oo) = (2/tt)Gf(u))) is given quantitatively by the portion arising from the 
motion of the solute and the single largest contributing solvent Moreover, as with the 
molecular dynamics studies, these findings seems to be completely independent of how oth- 
erwise collective the solvent motion is. The question we are left with, therefore, is whether 
the nature of the relaxation really does change qualitatively as we cross the band edge. It 
does not seem all that far fetched for this same two-particle scenario that works so well just 
inside the band to apply at higher frequencies as well. 



B. The High-Frequency Parts of the Vibrational Relaxation Correlation Function 



To make our subsequent development more concrete, suppose we specialize to the familiar 
case of a diatomic solute dissolved in an atomic liquid. Taking the atoms of the diatomic 
to be labeled as 1 and 2 (with masses m\ and m 2 ), and indexing the atoms of the solvent 
by J ' = 3, . . . , N, we find that if the potential energy of interaction between the solute and 
the solvent can be written in terms of site-site pair potentials 



v = M r ij) +M r 2j)} ; 
then the force acting on the diatomic's bond is given by 



(2.9) 



dV 

dx 



^12 X! 

i>2 



777.1 



m 2 



u 2 {r 2j )r 2j 



■ r i2 • 



(2.10) 



Here fi\ 2 = ni\m 2 j (m\ + 7772) is the solute reduced mass and r = r/|r|. The requisite time 
derivative of the bond force is thus 



F r 



i>2 

7>2 



— %(r y )ry(r lj ■ r 12 J 
mi 

— M i( r ij)4(rij -?i2) 
mi at 



1 

m 2 
1 

U 2^2j)-T + [j2j • ri2 J 

777 2 at 



u 2 [.r 2j )r 2j [r 2j ■ v l2 



(2-11) 



Let us now consider what the highest frequency parts of the force- velocity correlation 
function, Eq. ( |2.1| ), would look like for such a force. Generalizing our arguments in the 
previous section, we would expect that the dominant contribution would come (at most) from 
the two solvent atoms, call them a and b, nearest the two solute sites 1 and 2 (respectively) .0 
Furthermore, given the rapid variation of the site-site potentials ui^ 2 (r) in the region of 
interest, we might also expect the terms involving the second derivative u"(r) to be far more 
important than those involving the first derivative u'(r). Hence we should be able to express 
the correlation function as just 
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G F (t) w (^) (u"(r la (t))u"(r la (0)) cos9 la (t) cos0 lo (O) ha(t) fi„(0)) 
\mi/ 

+ f — V <u"M*)) w"(r 26 (0)) 008^(0 cos0 26 (O)r 26 (i)r 26 (O)) , (2.12) 

with cos 0Aj(t) = (f Aj -fi2) (^4 = 1,2). We could, in fact, make use of Eq. ( |2.12j) as it stands; 
since it requires that we know no more than the dynamics of 4 degrees of freedom, it already 
represents a major simplification over the full many-body problem. We can make one last 
simplification, however. The variation of the cos 6 Aj{t) angular factors is also guaranteed 
to be far smaller than that of the u"(r) factors at short times, if only because of the small 
dynamic range of the cosine function (—1 < cos# < 1). For the purposes of evaluating 
high-frequency behavior it therefore suffices to approximate the cosine by its t = value. 
Our final high-frequency expression for the vibrational relaxation correlation function thus 
reduces to a remarkably simple form 

G F (t) « (^) 2 (u"(r la (t))u"(r la (0)) cos 2 lo (O) r la (t) r la (0)) 

+ (^) 2 (u"(r 2b (t))u"(r 2b (0)) cos 2 6 2b (0)r 2b (t)r 2b (0)) , (2.13) 

relying only on the one-dimensional, purely radial dynamics of two solvent molecules moving 
with respect to the two solute sites. 



C. Mutual Nearest Neighbor Dynamics 

Of course, since we are in a liquid, the exact dynamics of even this small a number of 
degrees of freedom is still intimately connected with the dynamics of the remainder. Were we 
unable to achieve any reduction of this dynamics comparable to the simplification we found 
for the correlation function, there would thus be little practical advantage to our formalism. 
Fortunately, the same kinds of analysis we have been pursuing also reveal a wonderful 
dynamical separability: the important solvents end up being dynamically decoupled from 
the rest of the solvent as well. 

This basic notion of dynamical separability is illustrated by looking at the time evolution 
of a mutual-nearest-neighbor pair of atoms in an atomic liquid.ill£l The equations of motion 
for any two atoms a and 6, with masses m a and m&, are 



m„r, 



u'(r ab )r ab - u \ r aj)ij 



a L a " y abj'-ab / t >* y ajj'-ja 

m b r b = -u'(r ab )r ab - u \ r bj)^jb , 

provided all the atoms interact by pair potentials u(r). But, if these two particular atoms 
form a mutual-nearest-neighbor pair, then u'(r ab ) will be significantly larger than any of 
the other forces affecting the pair. Hence, to within a high level of accuracy, the relative 
acceleration will be given by 

r a b = n - r a w u {r ab )v ab , 
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meaning that the scalar displacement r ab will satisfy the differential equation 

r ab = —u'{r ab ) + r ab n 2 ab , (2.14) 



with fi = m a m b /(m a + m b ) the reduced mass of the pair. Here the angular velocity factor 

2 

ab 



VL 2 ah = — Tab ■ fab and we have used the vector identity 



?o6 • fab = r ab - r ab Q 2 



ab 



Equation (|2.14 ) is intriguing, but it is not as simple as it may appear; the evolution 



of the centrifugal acceleration term r ab VL 2 ab is controlled by other forces besides u'{r ab ) (and 
thus by the motion of other atoms besides a and b). However it is fairly easy to see that this 
term too is almost always going to be negligible compared to the intra-pair force. We can 
estimate its size by using the equipartion theorem to get the mean-square angular velocity^ 
VL 2 ah =< uj 2 >= [k B T I [ir 2 a} ). If we take the pair potential to be of the Lennard- Jones form, 
u(r) = 4e[(cr/r) 12 — (cr/r) 6 ], and we then make the (very) conservative estimate that a 
typical intra-pair distance r ab = a, we find that a typical value of the centrifugal force is 
(ksT l /e)(e/a), whereas a typical intermolecular force is 24(e/er). Hence under any liquid- 
state conditions, indeed at anything other than the most extreme supercritical conditions, 
we can safely neglect such centrifugal contributions.il The dynamics of a mutual-nearest- 
neighbor pair of atoms is therefore well described by a simple, purely radial, equation of 
motion: 

Tab = --u'(r ab ) , (2.15) 

even in a dense liquid. 

This same argument can easily be adapted to the vibrational relaxation example we 
considered in the last section. If only one atom (say number 1) of our diatomic solute forms 
a mutual nearest neighbor pair with a solvent atom (solvent a) , then the equations of motion 
for the critical solvent atom and for R cm , the position of the diatomic's center of mass, are 
given accurately byil 

m a r a = -u[(r la )r la (2.16) 

M 12 R cm = u[{r la )r la (2.17) 

where m a is the mass of the solvent atom, and M\ 2 = mi + rri2 is the total mass of the 
solute. For a rigid diatomic, though, 



so that the analogue of Eq. (|2.14j) becomes 



ha = — «'l( r la) + r la^L + ( TT~) r ^ha • ?i 2 , ( 2 - 18 ) 

A*a VMi 2 / 

with fi a = m a M 12 / (m a + M 12 ) the solute-solvent reduced mass. 
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As we did before, we can now turn to the non-radial terms. The second, solute-solvent 
centrifugal, term in Eq. (|2.18|) can be discarded for precisely the same reasons we just went 
through; on the average this term is far smaller than the radial force ^(r^). The third 
term seems to present a new case, inasmuch as it arises directly from the rotational motion 
of the diatomic solute. Yet, we can argue that it too should be small. Barring pathological 
situations, the interatomic separation and the reduced mass for the diatomic should be 
comparable to those of a solute-atom/solvent-atom mutual-near-neighbor pair. We should 
therefore expect similar results for the centrifugal terms 

^12 = ~ r 12 ' ^12 ~ = — f la ' ?la ? 

both of which are equal to angular velocity factors on the average 

'^12/ = ( f 12 • f 12 ) , (fl 2 la ) = (ha ■ ha 



and should therefore be small. By the same token, the average rotational term < — r la -f 12 > 
= < f la • f 12 >. But by Schwarz's inequality 



r la ' r 12 



< (f 12 -f I2 ) 1/2 (f.-M 1/2 



Hence, the radial force should dominate both non-radial terms, leaving us with the strikingly 
simple, one-dimensional equation of motion: 

ha = —u[(r la ) . (2.19) 

Ha 

In the event that both atoms of the solute are members of mutual-nearest-neighbor 
pairsE3 with solvent atoms (the situation described at the close of Sec. [11 B|, with the two 



solvent atoms labeled a and 6), similar manipulations and arguments tell us that the equation 
of motion should be generalized to the two coupled differential equations 

ha = — — w'i( r io) ~ i^—u' 2 {r 2b )cos6 ab 

Ha Wl 12 
1 1 

hb = u' 2 (r 2b ) - — w'i( r ia) cos9 ab , (2.20) 

Hb M12 

with fi b = m b M 12 / (m b + M 12 ). Nonetheless, consistent with our neglect of rotational dynam- 
ics, we can regard the relative orientation of the two pairs, cos 6^ = (f la • r 2b ), as a constant 
of the motion, so the resulting equations are still purely one-dimensional. We might also 
predict that, except for very rare liquid configurations, one of the two tagged solvent atoms 
will be far more important than the other in determining the force along the bond. We 
therefore expect that we should never have to do more than solve an equation of the form 
of Eq. ( |2.19| ) for the most important solvent/solute-site pair. 



D. Putting the Pieces Together 

Because all we need to know for a given liquid configuration is the location of the critical 
solvent closest to one end of the solute or the other, our working formula for the vibrational 
relaxation correlation function is simply, 
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G F (t) = - A p Jdrg mnn (r; ?ia) (u"[r(t)]u"[r(0)] cos 2 #(0) r(t)r(0)) r , (2-21) 

for any homonuclear diatomic dissolved in an atomic liquid. Here r is the distance between 
the special solvent and the nearest solute atom, <? mnn (r; f 12 ) is the probability density of r, 
what we shall call the solute's mutual-nearest-neighbor distribution function, and 8(0) is 
the initial orientation of the solvent/solute-atom pair relative to the bond direction, f 12 . 
The necessary dynamics launched from that configuration, r(t), is then prescribed by the 
equation of motion 

f(t) = —u'[r(t)} 

r(0) = r, f(0) = v, (2.22) 

with fi uv = mMj (m+M) the solvent-solute reduced mass (m being the solvent atom's mass 
and M being the solute's total mass). 

What liquid-structure effects there remain in the problem show up in the function 

9mnn{ r ] ^12Jj 

P9mnn(r; ?i 2 ) = ( E ' E ~ ~ *a)]) (2-23) 

\A=l,2 i>2, j=mnn of site A I 

which is the equilibrium probability density for finding a solvent atom that is both a mutual 
nearest neighbor and is displaced by r from one end of the solute (p being, the solvent 
number density). In the event that both solute sites are members of mnn pairs,t3 the primed 
summation selects the site, A, whose partner is closest. Consistent with this defintion, by 
the brackets with subscript r, we mean just the equilibrium average over initial velocities v 
in Eqs. Q2.21]) and ( 2.22|) , carried out at fixed r. 



The ultimate test of our approach, of course, will be its numerical accuracy in reproducing 
the true high-frequency vibrational friction, but even in advance of such tests it is worth 
looking at the internal consistency of the theory. Supposing for the moment that our basic 
idea is correct — that one or two critical solvents dominate the instantaneous friction — we 
still have to remember that the rest of the liquid continues to evolve while the critical solvent 
is moving. Is it even obvious that the particular solvent that makes up the initial mnn pair 
with the solute remains an mnn pair throughout the relevant time window? As we shall 
see, for Xe liquid the high-frequency part of the friction is largely determined by dynamics 
occuring at times less than 300 fs. A typical speed for the mnn distance is [k B T / ij, uv ) 1 / 2 , 
so at a temperature of 280 K, a diatomic such as I 2 dissolved in Xe could have its mnn 
distance change by about 0.5 A during this time. However, the mnn distribution itself has 
a full width at half maximum that is typically of this same order in dense XeB It therefore 
seems quite plausible that the identity of the mnn partner could be preserved long enough 
for the theory to operate successfully. 



III. INSTANTANEOUS NORMAL MODES AND INSTANTANEOUS PAIR 

THEORY 



The ease with which we can perform calculations with Eqs. (|2.21 ) and (|2.22 ) aside, 
we have not really finished our task; there is still a conceptual gulf we need to cross. Our 
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theory for high-frequency vibrational population relaxation, does not, in itself, resolve the 
central conflict between V-V and V-T pictures of the relaxation. On the one hand, the 
instantaneous-pair theory certainly seems to suggest that high-frequency vibrations are 
converted into the liquid's translational degrees of freedom. Indeed, the whole motivation 
for separating out the contributions of a single solvent is that the special solvent highest on 
the repulsive wall of the intermolecular potential completely overshadows any other solvent 
- and the fact that this sharply repulsive portion of intermolecular potentials is largely 
what comes into play in the evaluation of the equations of motion means that the dynam- 
ics is a case study in classical scattering rather than an example of bound-state motion. 
In short, the relevant dynamics seems to resemble translation more than intermolecular 
vibrations. By the same token, however, it is difficult to reconcile this stance with the 
noteable success of resonant-vibrational-energy-transfer models for frequencies within the 
liquid's band,0@ especially given the observations that there are binary features within 
the instantaneous-normal-mode band.0@ It would not be impossible for the mechanism of 
vibrational relaxation to switch as we go from one frequency regime to another, but such 
a pronounced dichotomy would be especially problematic in the immediate vicinity of the 
band edge. 

One way to think about this dilemma is to ask what would happen if we tried to con- 
struct an instantaneous-normal-mode treatment assuming that all of the assumptions of 
the instantaneous-pair theory were, in fact, correct. That is, suppose we assume that the 
vibrational-friction correlation really is given by the instantaneous-pair approximation of 
Eq. (|2.21[) and that the dynamics really ought to be that prescribed by the matching 



instantaneous-pair approximation of Eq. ( |2.22| ). Rather than solving Eq. ( |2.22| ) exactly, 
an INM theory would make an instantaneous harmonic estimate for the time evolution: the 
displacement from the initial pair separation would be written as0 



r(t) - r(0) ee q(t) = (— ^-A (1 -cos wot) + (— )smu t , (3.1) 



-^-2 J (1 - COS LU t) + ( — 

^u» w o / V No- 

where u ee \Ju"[r{Qi)]/ n uv is the instantaneous frequency of the "harmonic mode" and 
/o ee — u'[r(0)] is the instantaneous force. 

However, besides just assuming that the dynamics was harmonic, the traditional, lin- 
earized version of INM theory would go one, key, step further. The evolution of the coupling 
force, here u'(r), would be linearized in the mode coordinate^ 

u'[r{t)) w u'[r(0)]+it"[r(0)] q(t) . (3.2) 

This step is what leads to the appealing form of both INM linearized solvation theory and 
the INM theory for vibrational friction.&~il@ The appearance of what have been called the 
INM influence spectra stems from this very step. But it should be clear from our analysis of 
instantaneous-pair theory that this step is strongly suspect under the very conditions that 
lead to high-frequency relaxation. If the intermolecular potential is really sharply varying, 
then that is the last place where one should attempt a low-order expansion. 

We could raise the same kinds of doubts about the harmonic dynamics of Eq. ( |3.1| ). 
but suppose, for the moment, we pursue just the issue of nonlinearity in the coupling. How 
would the nonlinearity alone make itself felt? We can present this question more generally if 
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we return to a general solute-centered correlation function introduced in Sec. [II A| . Within 
the framework of the instantaneous-pair assumptions, the correlation function we need to 
evaluate is of the form 

G F (t) = (w'[r(t)\ w'[r(0)\ r{t) f(0)> , (3.3) 

where w(r) is some spectroscopic probe potential which may, but need not, be related to 
u(r), the pair potential governing the dynamics. 

Were we to linearize w(r) the same way that we did the force in Eq. fl3.2| ), Eq. ( p.l| ) 
would allow us to write the correlation function as the cosine transform of an influence 
spectrum, p(ou), a standard INM outcome. 

Qhnear-INM^ = J duJp ^ COSUJ t. (3.4) 

PH = — (K) 2 5^-^)) • (3.5) 



Here [i uv is the appropriate reduced mass, w' = u/[r(0)], ujq = yu"[r(0)]//j, uv , and the 
brackets denote the (fully anharmonic) equilibrium average over the initial coordinate r(0). 
As we might have expected, however, this influence spectrum is too limited; it spans a range 
of frequencies no larger than the INM band (the density of states) itself 

D(u) = (6(u-w )) . (3.6) 

We want to point out, though, that the instantaneous-pair formulation of the problem makes 
it simple enough to include the full nonlinearity of the coupling, even while retaining the 
harmonicity of our dynamics. If we carry out the equivalent of Eq. ( |3.2| ) to all orders 



00 1 

= £ q n (t) , (3.7) 



n=0 U - 



with Wq 1 ^ the n-th derivative of w (r) evaluated at the instantaneous position r(0), we can 
still substitute Eq. (|3.1|) into Eq. (|3.3|), providing us with an expression for the necessary 
correlation function in terms of an infinite sum 



G™(i) = J2 G F(t) , (3.8) 

n=l 

C$\t) = (4^ <f-\t) q(t) q(0)) , (3-9) 

though the bracket now represents the average over the initial velocity as well as the anhar- 
monic average over the initial coordinate. But, as we show in the Appendix, this velocity 
average can be performed analytically, yielding the result that the n-th term is 

GP(t) = ^ (W V n [t- r(0)]) , 

H"uv 
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V n [t; r(0)] 



•/Vi7) n 



a 



{(ncoauot)H n ^[iX(t)] + (n - l)[2iX(t)]H n _ 2 [iX(t)]} 



X(t) = ^— j tan— , Y(t) = —. (3.10) 

In these formulas H n (x) is the n-th order Hermite polynomial^ and we have defined the 
chararacteristic length scale and time scale for the dynamics by the instantaneous-position- 
dependent quantities 

a " = — ' t0= [k^i) • (311) 

respectively. 

The first term in this series is, of course, the standard linearized INM result, 

G$\t) = G h F near - INM (t) , (3.12) 

involving a configurational average of (something proportional to) coswoi, but each higher 
order term clearly brings in higher powers of the trigonometric functions. In fact, since 
H n (x) = (2x) n plus lower powers in xji it is easy to show that each V n integral is a polyno- 
mial form in coso^ and sinwot, and that the leading (highest order) term is proportional 
to cos n ujpt = cosnuJot + lower order harmonics of ujq. In other words, each successive term 
brings in successively higher-order overtones of the fundamental INM band of the liquid. 
The instantaneous pair theory with INM dynamics is thus a genuine multi-phonon theory 
of high-frequency vibrational relaxation in liquids — as long as the proper respect is paid 
to the nonlinearity of the coupling.E3'E3 

We should hasten to point out that a literal term by term evaluation of this sum is 
not necessarily going to be of any great computational help. For one thing, it is not 
straightforward to write the n-th order equivalent of the (one-phonon) influence spec- 
trum, Eq. (|3.5|). Since V n not only has terms of the form cosnwot, but can include 
cos (n — l)uot, cos (n — 2)u t, . . . , and cosu t terms as well, the final n-phonon spectrum 
could have contributions from V n , V n+ i, and all of the higher-order integrals. A more fun- 
damental way of making this same point is to remember that there is no real reason to 
expect that any low-order truncation of the multi-phonon series will adequately describe 
the relaxation process. A more realistic approach to the problem might require us to sum 
the series to infinite order, at least in some approximate fashion. 

It is straightforward enough to evaluate the combination of Eqs. fl3.1| ) and fl3~3] ) nu- 
merically, without making any approximations, but it is revealing to examine a special case 
for which we can perform the necessary multiphonon sum analytically. If the spectroscopic 
probe potential is an exponential repulsion (a fairly common assumption in theories of vi- 
brational relaxation)!ll!el 

w[r(t)] = w[r(0)] e~ a[r{t) - rm , (3.13) 



with a a constant, then all of the derivatives in Eq. (|3.10|) can be expressed in terms of a 



and w' = w 



(i) 
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G F (t) = (—)((w>) 2 V[t; r(0)]) 



V[t; r(0)] = r(0)] K —^— . (3.14) 

72=1 \ n l )- 

As detailed in the Appendix, the sum can now be evaluated in closed form in terms of the 
ratio 70 = a/a . 

V[t- r(0)] = S[t; r(0)] e R[t ' r(0)] , 

r(0)] = 7o(cos^ t-l) + 7o 2 (^r~) , 

5*[t; r(0)] = cost^ot — 7osin 2 u; t + 27 2 cosu^ f ~ — _ (3.15) 

In fact, for vibrational relaxation the situation becomes even simpler.il If the interatomic 
potential itself is exponential 

u[r(t)] = u[r(0)] e' a[r{t) - rm , (3.16) 

then w(r) = u'(r) is proportional to the same exponential. The ratio 70 = a/a^ is thus 
identically equal to 1 and the charateristic time reduces to the constant 



Equations ( |3.14|) , ( |3.15| ), and ([3.17]) constitute our basic non-linear INM theory for 



high-frequency vibrational relaxation. We shall see presently how well it fares compared 
to nonlinear-coupling plus anharmonic-dynamics treatment of the full instantaneous-pair 
theory. 

IV. NUMERICAL RESULTS 

A. Numerical Methods and Models 

In this section we numerically evaluate the vibrational friction felt by a rigid homonuclear 
diatomic molecule dissolved in an atomic liquid. The atoms of the solvent are taken to 
interact with each other and with the two atoms of the solute through Lennard- Jones pair 
potentials 



u(r) = 4e 



r 



though the solvent atom-solvent atom (yv) and solute atom-solvent atom (uv) Lennard- 
Jones parameters can differ, as can the masses of the solute and solvent atoms (m u and m v , 
respectively). While the effect of the solvent on the intramolecular vibration is the whole 
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point of these calculations, for the purposes of computing the friction, the solute atoms can 
be regarded as being separated by a fixed distance r eq . 

For purposes of defmiteness, we regard our solvent as fluid Xe (a vv = 4.1 A, e vv = 
222 K, m v = 131.3 amu).^l The two specific solution models we consider are the standard 
Tuckerman and Berne model,0'0 in which the two solute atoms are regarded as identical to 
the solvent atoms, 

and the model used by Brown, Harris, and Tully,0 

a uv = 3.94 A, e uv = 324 K, m u = 126.9 amu, r eq = 2.67 A . 

This latter model was explored largely in response to experimental studies by Paige, Russell 
and Harris of the relaxation of highly vibrationally-excited I2 dissolved in Xe.0 

Our theoretical predictions were constructed so as to rely solely on the equilibrium prop- 
erties of the solutions being studied. We were therefore able to obtain the information we 
needed by using standard canonical-ensemble Monte Carlo simulations^ with 106 solvent 
atoms and an immobile diatomic solute. All of these simulations employed periodic bound- 
ary conditions and had their maximum-move distances adjusted to achieve a 50% acceptance 
rate for the Metropolis attempted moves. Ensemble averages were then computed by sam- 
pling fluid configurations every N t such attempts, with N t chosen so as to ensure minimal 
configuration-configuration correlation in the force along the diatomic's bond. The actual 
values of N t ranged from a low of 2000 for dense gases to as high as 5000 for high-density 
fluids. When exact dynamical results were needed to give us something with which we could 
compare, we carried out full molecular dynamics simulations on the same rigid-diatomic- 
plus- 10 6-sol vent system, now allowing the diatomic to translate and rotate. Specifically, 
classical trajectories of length tlj [tlj = (m vv a^ v / 'e^,) 1 / 2 = 3. A7ps] were computed start- 
ing from each of 10 5 (uncorrelated) Monte-Carlo generated configurations using the Rattle 
version of the velocity-Verlet integratorEl with a time step of 0.001 tlj. The initial velocities 
for these trajectories were selected from the Maxwell-Boltzmann distribution.0 

Results for the pair theories were obtained in much the same way as the exact molec- 
ular dynamics, but the requisite trajectories involved only the motion of one-dimensional 
solute-solvent distances r(t). The outcome for the full instantaneous-pair theory could 
therefore calculated by using a simple velocity-Verlet integration starting from the Monte- 
Carlo generated liquid-state configurations.^" Averaging over the initial velocities was then 
accomplished by using 40 velocities equally spaced between ±(40A;bT/ / u u1 ,) 1 / 2 , appropriately 
weighting by the Maxwell-Boltzmann distribution. The computation with the alternative 
pair theory, the corresponding nonlinear-INM approach, was even easier, since it only needed 
the analytical dynamics of Eq. ( |3.1| ), though it too required numerical averaging over the 
initial velocities. 

All cosine transforms were performed using standard FFT methods.El 



B. Basic Results 

We begin our exploration by looking at the vibrational friction felt in the Tuckerman and 
Berne model of a diatomic solute dissolved in a high density supercritical fluid. As one can 
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see from the bottom panel of Fig. |I|, our earlier linearized-INM theory (the long-dashed line) 
does a reasonably credible job of reproducing the shape of the exact-molecular-dynamics- 
derived friction (the solid line) within the band of the liquid. However once we venture 
beyond the band edge (upper panel) the discrepancy between the two becomes painfully 
obvious. By contrast, the full instantaneous pair theory (the short-dashed line) provides 
a virtually quantitative match to the exact results at these higher frequencies. Apparently 
combining a nonlinear treatment of the solvent coupling with an anharmonic approach to 
the solvent dynamics amply suffices to capture the missing high-frequency response of the 
solvent — and it does so despite our omission of all but an infinitesmal fraction of the solvent 
dynamics. The relatively poor performance of the same instantaneous-pair theory within 
the solvent band (lower panel) continues to remind us, at the same time, that an accurate 
treatment of the low-frequency regime requires a proper inclusion of the more collective 
aspects of the solvent's motion. 

These same points are made more explicitly by using Eq. ( |1.1| ) to compute the actual 
vibrational energy relaxation rates from the friction. As we can see from Table pj the 
instantaneous-pair theory precisely mimics the two-orders-of-magnitude drop in relaxation 
rate that transpires once we go from vibrations in the middle of the liquid's band to those 
100 cm -1 outside the band. Yet, within the band, the rate predicted by this theory is 
consistently too small, just what we might have surmised would happen if we were leaving 
out the dynamical channels for relaxation provided by many-body motions. Linearized INM 
theory, on the other hand, behaves in a perfectly complementary fashion. It gives sensible 
relaxation rates within the band, but can do no better than trace out the precipitous decay 
of the band edge once we start to leave the band. 

Our next question, clearly, is whether the success of the full pair theory means that a 
harmonic approach to the solvent dynamics suddenly becomes inappropriate at high fre- 
quencies. Looking again at the high-frequency friction, Fig. ||, we find that if retain the 
exact nonlinearity in Eq. ( |2.21| ), but use instantaneous-normal-modes to prescribe the 



mutual-nearest-neighbor-pair dynamics, Eq. ( |3.1| ), the resulting friction spectrum is still 
virtually quantitative in its agreement with the exact molecular dynamics. As we might 
have anticipated from the multiphonon discussion in the last Section, having nonlinear cou- 
pling allows even harmonic dynamics to span a realistic frequency range. What we might 
not have predicted, however, was that the coupling was so much more important than the 
anharmonicity that incorporating the nonlinearity of the coupling was all one needed to do 
to recover accurate vibrational relaxation rates. 

The level of similarity between these two pair theories is actually even greater than it 
might appear on the scale of Fig. One of the more dramatic findings of the studies of 



vibrational relaxation in the solid state is that once the amount of e 



rergy transferred gets 



large enough, relaxation rates tend to obey an exponential gap law.iJliio That is, the rates 
depend on the frequency of the solute vibration uj roughly as 

Tf 1 ~ exp-(^— ) , (w» u char ) , (4.1) 

where uo c h ar is some characteristic phonon frequency of the solid. There is a long history 
to these kinds of exponential relations for a variety of different kinds of energy transfer 
processes,^ but the particular application to solids is well known to follow directly from 
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a multiphonon picture of the relaxation: steepest-descent evaluation of the relaxation rate 
predicts the rate to be exponential in the number of phonons excited. 

What happens in liquids?0S A logarithmic plot of the relaxation rate, Fig. |3], shows that 
the decline in the molecular-dynamics-computed rate with solute frequency is indeed close 
to exponential, a trend reproduced exceptionally well by the instantaneous pair theory. In 
fact, given the numerical difficulties inherent in extracting the high-frequency components 
of a molecular-dynamics trajectory^ we might venture to suggest that the pair theory 
results are actually the more reliable of the two at high frequencies. 0@ Interestingly, the 
exponential line that the two curves seem to follow is 

Tf 1 ~ exp-(^o) (4.2) 



with to the very time scale prescribed by nonlinear INM theory, Eq. ( [3.11| ). This statement 
actually requires a little qualification. Our formula for t expresses it as a function of cto, 
which, itself is a function of the initial mnn distance r(0). However we can evaluate an 
average «o a priori by using the «o values of Eq. fl3.11|) for the solute-solvent Lennard- 
Jones potential averaged over the equilibrium distribution of initial mnn distances, Eq. 
( |2.23| ). The resulting a is what is used to generate the t for the long-dashed line in the 
figure. El 

Of course, coming back to the solid-state analogy, we could regard an exponential gap law 
coming out of nonlinear INM predictions as completely consistent with our multiphonon-like 
interpretation of the INM theory. In fact, it is not hard to show that our simple analytical 
results of Eqs. (|3.14j) , ( |3.15| ), and (|3.17| ) do predict the same exponential decay as the full 
instantaneous pair theory (Fig. |3|, bottom panel). Apparently, even with this intrinsically 
vibrational picture of the solvent dynamics embedded in the calculation, the agreement 
between the pair theories with harmonic and anharmonic dynamics is robust enough to 
span a number of decades of solvent response. 

It is revealing to think in some detail about this remarkable agreement between the 
full instantaneous pair theory and its nonlinear-coupling/harmonic-dynamics analogue. A 
closer inspection of the two theories in the time domain, Fig. [|, emphasizes the apparently 
fundamental differences between the two views of the dynamics. If we restrict ourselves 
to a single initial liquid configuration, as the figure does, we see that the harmonic curve 
reflects the oscillatory character of bound-state motion, whereas the lack of periodicity 
in the full (anharmonic) result is consistent with what the pair dynamics really is — a 
scattering process. So why should such disparate views lead to the same results? The 
answer is suggested by the third curve in Fig. 4, which also assumes harmonic dynamics but 
truncates the coupling at linear order, precisely the sort of theory we would configurationally 
average to arrive at Eqs. ( |3.4| ) and Q3.5| ). It is evident that neither of the two harmonic 
curves provides all that faithful a rendering of the anharmonic dynamics, but what stands 
out is that the nonlinear coupling allows for a much more accurate reconstruction of the 
correlation function at short times, reproducing not only the initial decay of the exact 
correlation function, but also its one oscillation.^ 

It certainly makes sense that the highest frequency behavior of the friction depends 
critically on the shortest time dynamics, but we can say a little more than that. Obviously, 
to get the correct high-frequency Fourier transform of Gp{t), it is important that we know 
its behavior everywhere it changes rapidly — which includes both the initial fall-off and 
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the subsequent rebound we just referred to. Equally important, though, is the fact that 
when these single-configuration correlation functions are averaged, the spurious oscillatory 
behavior of the harmonic picture that we see at longer times, which seems to portray solvent 
molecules as colliding with the solute over and over again, dephases away to nothing. Thus, 
at least within the time window provided by the rapid decay of the mnn bond-force- velocity 
correlation, vibrational motion can indeed manage to look identical to single scattering 
events. To put it another way, for our own limited purposes of computing vibrational 
population relaxation rates, we can evidently continue to regard the essential liquid dynamics 
as intermolecular vibrations, both at high and at low frequencies. 



C. Thermodynamic State Dependence and the Application to I2 in Xe 

We close our discussion with a look at a more experimentally oriented calculation, that for 
the vibrational population relaxation of the first vibrationally-excited state of I2 dissolved 
in fluid Xe at 280 K.0H 

To begin with, we should note that from xenon's perspective, the I2 vibrational fre- 
quency of (oo/2ttc) = 211 cvrC x actually qualifies as a rather high frequency. As is clear 
from Fig. |5], this frequency lies well outside the INM band of Xe liquid (and even further 
outside the band of lower density Xe). It therefore confronts us with precisely the issues we 
have been discussing. Another interesting feature of this example, though, is that at this 
temperature, Xe can be taken from dense liquid, through the liquid-gas coexistence region, 
into a rarified gas. We can therefore examine a range of dynamical possibilities within this 
one example. The only additional step we need to take is to provide some way to compute 
the relaxation rates within the coexistence region using data taken only from single-phase, 
constant temperature and volume, simulations. 

To meet this last requirement, we assume that the rate at which a solute relaxes vibra- 
tionally, T _1 , is far faster than the rate at which it moves between liquid and gaseous regions. 
Under these conditions, the experimentally measured relaxation rate should be simply the 
weighted sum 



+ * 9 as(T-i) gas , (4.3) 



where {T~ l )u q and (T _1 ) gas are the relaxation rates calculated along the coexistence line 
in the liquid and gaseous phases, respectively, and xu q and x gas are the respective mole 
fractions of the two phases at the given temperature and total number density p. These 
mole fractions, in turn, are derived from the lever rule 

•Eliq 1 Xgas > ( / ^ -Z ^) 

Pliq Pgas 

and the values of pu q and p gas , the individual densities of the liquid and gas phases. 

Accordingly, if we use the full instantaneous pair theory to evaluate the relaxation rates 
in the liquid and the gas, and we allow Eqs. ( fOl) and (|4.4|) to interpolate between the 
liquid (1.7 g/mL) and gaseous (0.44 g/mL) coexistence densities,^ we find that we can 
follow the vibrational lifetime of I2 over a wide range in solvent density (Fig. ^|). The 
general behavior shown here, that vibrational lifetimes decrease monotonically with density, 
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is hardly surprising, but a number of features of this graph are worthy of some attention. For 
one thing, we note that the simple pair theory predicts lifetimes that are essentially identical 
to those computed by the elaborate molecular dynamics simulations of Brown, Harris, and 
Tully, at least at the two densities that those authors were able to study.0 However, the 
computational simplicity of the pair theory makes it straighforward for us to trace the 
overall trends with thermodynamic state as well. The basic increase in relaxation rate that 
we see with density (a 16-fold increase over the same interval where the density increases 
by roughly a factor of 6) is something that one would have expected based in independent- 
binary-collision (IBC) models.il'ill! Those models, though, would attribute the increase to 
a hypothetical growth in the "collision rate". In the instantaneous-pair theory, by contrast, 
the solvent density influences only the equilibrium probability of having a given mutual- 
nearest-neighbor pair distance; higher densities translating into smaller average distances. 
The subsequent dynamics — and therefore the relaxation rate for that mnn distance - 
is controlled solely by the temperature, which sets the initial velocity distribution of the 
solute-solvent pair, but which we hold fixed throughout the calculation here. 

The similar experimental consequences of the two different theoretical models make it 
difficult to distinguish them just by measuring relaxation rates versus density. We empha- 
size, however, the IBC model's reliance on an ill-defined collision-radius parameter,@I!Il 
something the more microscopically defined instantaneous-pair theory has no need for. In 
all fairness, we should also point out that both models lead to results that are extraordinarily 
sensitive to the precise shape of the repulsive wall of the solvent-solute potential. Indeed, the 
potentials employed here (and by Brown et al)0 predict relaxation rates differing by about 
a factor of 5 from the experimental estimates.0 That the likely origin of this discrepancy 
does lie in the intermolecular potentials is supported by our finding that small changes in 
the repulsive portion of the potential can shift our relaxation times by as much as an order 
of magnitude. This sensitivity implies that we could tune our results to be more in accord 
with the experimental data,t3 but since such curve fitting is not particularly germane to the 
issues at hand, we will leave the matter here. 



V. CONCLUDING REMARKS 

The overriding goal of this paper was not really to find out how quickly dissolved 
molecules can get rid of large quantities of vibrational energy, but to learn how we should 
visualize such processes in molecular detail. To that end we deliberately phrased the issues 
in terms of somewhat oversimplified caricatures. Genuine molecular motion in liquids is nei- 
ther the simple vibrations of crystals nor the occasionally interrupted ballistic motion found 
in gases. What we envisioned, though, was that keeping these two limiting cases in mind 
would serve to frame the possibilities. Even if the literal truth ended up being somewhere in 
the middle, we presumed that we would still be able to say whether molecules tended to lose 
vibrational energy by finding a matching vibration-like motion somewhere in the solvent 
(a more-or-less V-V event) or whether it was better to regard relaxing solutes as simply 
adding generic translational kinetic energy to the solvent (mostly a V-T process). 

Much to our surprise, the real situation turned out to be far more narrowly defined than 
we anticipated. In effect, we found that the structural disorder of liquids conspires with 
the short-ranged, rapidly-varying character of repulsive intermolecular forces to prescribe 
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a remarkably specific route for vibrational relaxation. The key to our understanding of this 
specificity was the idea of taking an instantaneous perspective: of looking at each liquid 
configuration and asking what the subsequent time evolution ought to be for short time 
intervals. The sharply varying character of the largest magnitude intermolecular forces tells 
us that out of all of the solvent molecules that could participate in the relaxation at a given 
instant, only the very nearest solvent should play a significant role in determining the high- 
frequency component of the solute-solvent coupling. This observation is central because it 
means that we only need to ascertain the dynamics of the solute and the special solvent. 
What is both more important and more unexpected though, are two realizations stemming 
from the basic properties of liquids: first, when the liquid structure is such as to have the 
special solvent molecule closer to the solute than to any other molecule in the solution (i.e., 
when the solute and the special solvent molecule are a mutual-nearest-neighbor-pair), then 
the pair motions themselves will dynamically decouple from those of the rest of the system; 
second, that these particular liquid-disorder-induced situations are precisely the critical 
ones dominating the high-frequency relaxation dynamics. For the kinds of intermolecular 
distances and potentials typically found in liquids, then, we see that the crucial dynamics is 
simply that of a two-body scattering process. 

Basic to all of this discussion is our assumption that the critical portion of the time 
evolution takes place at short times. In part, this assumption is valid because of our concern 
about high frequencies,!^ but equally, it stems from the short (few hundred fs) decay times 
characteristic of the force-autocorrelation function which serves to determine the vibrational 
friction. Of course, this time scale, as well, has its origin in the same fundamental ideas 
that we have been discussing. The harsh repulsive forces lead to a rapid decorrelation of the 
molecular trajectories and the average over the structural disorder of the liquid washes out 
any longer-time correlations that might remain. 

The presence of this short time scale was important to us not only in that it let us 
formulate our theory, but in its implications for the choice between V-V and V-T mecha- 
nisms. The fact that the dynamics has a scattering flavor is certainly more reminiscent of 
the gas phase than the solid state, but what we discovered was that for times this short, 
the distinctions simply evaporate. There is no discernable difference between the high- 
frequency vibrational friction evaluated as a scattering problem and the friction evaluated 
with instantaneous-harmonic — vibration-like — bound-state motion. As long as the cou- 
pling is sufficiently nonlinear, harmonics of the fundamental vibration can evidently play 
the precise role occupied by more impulsive kinematics. The short answer to our basic 
question is therefore that we can always think of solute as resonantly transfering energy to 
vibrational motions of the solvent; to the instantaneous normal modes for solute energies 
within the band of the liquid and to overtones of those modes for higher energies. This last 
interpretation is obviously not unique, but it has the conceptual advantage of allowing us 
to think of the progression from lower to higher vibrational energies as a gradual evolution 
in the kinds of dynamics the liquid has to bring to bear on the relaxation. 

This continued success of harmonic perspectives on liquid dynamics actually sends yet 
another message by highlighting one of the issues we raised in the Introduction. It is hardly 
disputable that understanding anharmonicity is vital to understanding liquid dynamics, but 
the distinctions between the different flavors of anharmoncity have significant consequences. 
If we needed to be able to model the detailed breakdown of the harmonic-mode picture of the 
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underlying dynamics in order to treat vibrational relaxation, we would have had a daunting 
task in front of us. Indeed, before we launched this work, this worrisome prospect was a 
legitimate concern; over the long times that it actually takes a molecule to rid itself of its 
vibrational energy, individual instantaneous normal modes completely lose any semblance 
of their orig inal MentitiesB What we have discovered, though, is that the primary issue 
in high-frequency behavior is the nonlinearity of the coupling being driven by the liquid 
dynamics, not the nonlinearity of the dynamics per se — and that this much more limited 
variety of anharmonicity is not all that difficult to understand. 

This same kind of distinction is actually a rather familiar one to spectroscopists. Early on 
in the history of infrared and Raman spectra, distinctions were drawn between "mechanical" 
and "electrical" anharmonicities in the vibrational spectra of individual molecules.il More 
recently, however, the issue has been the focus of considerable attention in the context of 
trying to understanding fifth-order nonresonant Raman spectra of liquids as a wholeS'ii^ii 
Purely harmonic systems would give rigorously zero signal in these experiments, so the very 
existence of a measurable signal is evidence of anharmonicity. But as Okumura and Tan- 
imura were quick to point out ,00 the same signals could arise equally well from fundamental 
liquid-dynamical anharmonicities or from nonlinearities in the polarizability dependence on 
normal modes (violations of the Placzek approximation) .il Just as we found here, the con- 
tribution of derivatives beyond the first for a spectroscopic probe potential (whether it be 
a force on a bond or a many-body polarizability) can be quite effective in extending the 
ways a liquid can respond to a probe. In fact, in perfect analogy to what we saw here, 
the polarizability nonlinearities can apparently be a much larger determinant of the spectra 
than the underlying liquid anharmonicities.^! 

Curiously, the results of this paper seem as if they can also be taken as supporting 
the independent-binary-collision model of vibrational relaxation. The centrality of solute- 
solvent pairs in the process was, if anything, put on firmer microscopic ground by our results. 
However, there are some key distinctions between IBC theory and what we have presented 
here. The most important of these is where the equilibrium considerations stop and where the 
dynamical concerns take over. Within IBC theory the existence, the precise definition, and, a 
fortiori, the rate of collisions all rely heavily on the liquid-state setting. Only the collision- 
induced vibrational relaxation itself is thought to be independent of the medium.0ii3~H 
With our instantaneous-pair model, though, the sole role of the many-body environment is 
to prescribe the equilibrium distribution of special (mutual-nearest-neighbor) solute-solvent 
pairs. These pair distances form the initial configurations for the dynamical calculation, but 
all of this subsequent dynamics, including the distribution of initial pair velocities, are purely 
few-body in character. We should emphasize that these distinctions are not merely technical. 
The arbitrariness of IBC theory results from the difficult task it has of solving enough of 
the nonequilibrium statistical mechanics of the solution to extract a cleanly separable pair 
motion from the remaining background. Lacking a definitive solution, IBC applications are 
compelled to adopt ad hoc definitions of collisions, postulating such criteria as minimum 
collision radii.il'il By contrast, the crisp separation with instantaneous pair ideas between 
many-body equilibrium information and few-body dynamics, not only removes any need 
for guesswork, it might even allow for future quantum mechanical treatments of vibrational 
relaxation based on the relative ease with which we can do two-body quantum scattering 
calculations.0 
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There are some other directions the instantaneous pair theory should be extended in 
as well. All of the numerical examples explored here were based on the rather limited 
possibilities presented by a diatomic solute dissolved in an atomic solvent. The easiest 
generalization of this work to polyatomic solutes and solvents would probably have us focus 
on mutual nearest neighbor pairs between individual solute sites and solvent sites. However 
the success of such an approach is far from assured; sites within a molecule are always 
going to be strongly correlated. Thus in spite of the strongly binary character already 
seen in INM vibrational friction spectra with molecular solvents,!! a fully general version 
of instantaneous-pair theory remains to be demonstrated. The potential applicability to 
understanding how solvents mediate the vibrational energy relaxation of complex multi- 
mode polyatomics, in particular, presents a natural extension we find especially intriguing. 
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APPENDIX A: ANALYTICAL CALCULATIONS WITH NONLINEAR INM 

THEORY 



The velocity average needed in Eq. ( |3.9| ) is a Boltzmann average over the initial velocities 
v = q(0) 



V n (t) = 



q n -\t) q(t) q(0) 



(Al) 



with the mode displacements q(t) specified by Eq. ( |3.1| ). The value of such integrals will 
depend on the initial coordinate r(0) as well as on time, (as indicated by the notation in 
the text), but we have suppressed this dependence here for notational simplicity For the 
purposes of this appendix, brackets will always refer to averages over initial velocity. 

Evaluating expressions of this form is performed most easily by observing that we can 
calculate the average of the mode displacement raised to any power 



([oro q(tT) = [-i Y (tT H m[i X(t)] 
in terms of the Hermite polynomials H m (x)j^ 

H (x) = 1, #i(x) = 2x, H 2 {x) = Ax 2 



(A2) 



2,... 



the constants a and t defined in Eq. ( |3.11| ), and the functions X(t) and Yit) specified in 
Eq. ( |3.10|) . Equation (|A2|) results quite naturally from the integral relations 



7T 



I / dx (x + z) m e~ {ax) 



(2ia) 



H m (iaz) 



(A3) 



To have Eq. ( |A2| ) help us with Eq. (^T|), though, we need to express the remaining 
factor in the average in terms of the "INM basis", that is, in powers of q{t). After some 
algebra, we find 



q(t) q(0) 
Z(t) 



' k B T \ 

y (J"uv / 
COS Ld t 



X 2 (t) {Z(t)[a q(t)} 2 + [a q(t)] - l} 



^1 — COSCUot) 2 



(A4) 



Substituting Eq. ( |A4j) into Eq. (|A1|), repeatedly making use of Eq. (^), and employing 
the Hermite polynomial recursion relation!^ 



H, 



n+l 



2x H n (x) — 2n H n _Ax) 



then leads us directly to the expression for V n (t) given in Eq. (|3.10| ). 

To evaluate the multiphonon sum we need to sum these integrals to all orders — which 
we can actually do for the exponential model, Eq. Q3.13j ). For this particular model, the 
first two contributions to the multiphonon sum V(t) given in Eq. (|3.14j ) are 



-a) 

"of 



Vi(t) 



cos u t 
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-*y 

1! 



V%(t) = 7o(cos2a;ot — cosu t) , 



which makes for a particularly simple form for the linear and linear-plus-quadratic order 
terms in the vibrational relaxation case (70 = 1) 



V(t) linear = cos cu t 

y ^linear + quadratic = cos2uQ t ; 

results that clearly emphasize the one- and two-phonon character of these contributions. 
Higher order term do not preserve this special structure, but we can still perform the full 
sum by taking advantage of the generating function for Hermite polynomials:!^! 



00 , n 



n=0 



TV. 



HJx) 



exp (2xz 



Application of this formula and a simple variant lead to the final result given in Eq. 15| ) . 
The same result could also be obtained much more directly by simply perfoming the velocity 
average in Eq. ( |3.3|) . 
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TABLES 



TABLE I. Rates of Vibrational Population Relaxation for the Tuckerman-Berne Model System' 



lu/2ttc (cm- 1 ) b 1/Ti (ps^ 1 ) 

IP Theory 6 MD d linear INM Theory 6 

31 L3 O L2 

46 1.2 5.7 6.6 

77 0.85 2.9 4.7 

154 0.23 0.35 3.0 xl0~ 2 

231 4.5 xl0~ 2 4.4 xl0~ 2 2.2 xl0~ 7 

270 1.7 xl0~ 2 1.5 xl0~ 2 / 



a Landau-Teller-theory calculations for a vibrating diatomic molecule dissolved in a super- 
critical atomic liquid at reduced density pa 3 = 1.05 and reduced temperature /cgT/e = 2.5. 
Numerical values are those appropriate for a solute in a Xe solvent; to convert rates to those 
for an Ar solvent, both oo and 1/Ti need to be multiplied by Tx e / T Ar — 3A7ps/2.16ps = 1.61. 
b Vibrational frequency of the isolated diatomic solute. 
c Predictions from the full instantaneous pair theory. 
d Exact (Landau-Teller) calculations from molecular dynamics. 
e Predictions from linearized instantaneous-normal-mode theory. 

f This rate is zero to within the numerical accuracy with which we can determine the mag- 
nitude of the INM density of states. 
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FIG. 1. The vibrational friction felt by a vibrating diatomic molecule dissolved in a super- 
critical atomic liquid. Shown here are a number of different calculations for the real part of the 
frequency-domain vibrational friction of the Tuckerman-Berne model system (under the standard 
model conditions of reduced density pa^ v = 1.05 and reduced temperature fe^T '/e vv = 2.5) in units 
appropriate to a Xe solvent. The lower panel is the friction relevant to diatomic whose vibrational 
frequency falls within the expected bandwidth of the solvent; the upper panel portrays the friction 
for frequencies above the band edge. (Note the rather different scales used in the two panels.) 
In both panels, the rough solid line results from an exact molecular dynamics calculation of the 
bond-force autocorrelation function, the long dashed line (visible only in the lower left hand corner 
in the upper panel) shows the analogous predictions from linearized instantaneous-normal-mode 
theory, and the short dashed line gives the analogous predictions from full instantaneous-pair 
theory. 
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FIG. 2. The vibrational friction felt by a high-frequency vibrating diatomic molecule dissolved 
in a supercritical atomic liquid. As in the upper panel of Fig. [I], the calculations are performed for 
the standard Tuckerman-Berne model system (under the same thermodynamic conditions) with Xe 
as the solvent, and the solid and short-dashed lines show the full-molecular-dynamics and the full 
instantaneous-pair theory results, respectively. Here, however, the long-dashed line represents the 
outcome of the nonlinearly-coupled instantaneous-normal-mode version of the pair theory. Both 
pair theories interpret the relaxation as a few-body event in which only the initial conditions are 
prescribed by the fluid. They differ in that the full instantaneous-pair-theory uses the anharmonic 
dynamics of the key atoms to drive the evolution of the nonlinear solute-solvent coupling, whereas 
the INM equivalent assumes harmonic dynamics drives this coupling. In contrast to both of these 
treatments, the linear INM theory of Fig. [l] limits itself to linear coupling. 



36 



O 

X 



10 



10 



10 



10 



10 



10 



10 



10 



10 





150 



225 



300 



375 



450 



(o/2nc (cm ) 



FIG. 3. The vibrational friction felt by a high-frequency vibrating diatomic molecule dissolved 
in supercritical atomic Xe. In the upper panel, the instantaneous-pair (short dashes) and molec- 
ular dynamics (solid line) results from Fig. ||| are replotted here on a semilog scale to emphasize 
the near exponential-gap-law behavior. The lower panel compares the instantaneous-pair theory 
of the upper panel to the analytical results, Eqs. ( 3.14[ ), ( 3.15| ), and ( 3.17 ), where to generate 
to we have used the value of ao found by averaging over initial mnn distances. For the super- 
critical fluid, this average ao is 17.2 a^ e = 4.2 A" 1 , which yields a characteristic timescale, 
to = 0.042 tlj = 146/s. Numerical noise clearly makes it difficult to follow the exact molecular 
dynamics results over more than a limited frequency interval. It is only by going beyond this 
interval, however, that the small deviations from strictly exponential behavior begin to manifest 
themselves in the form of a slight downward curvature to the predictions of the pair theories. 
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FIG. 4. The bond-force-velocity autocorrelation function for the standard Tuckerman-Berne 
model system evaluated at a single initial liquid configuration for three different pair theories. 
Shown here is a normalized version of the instantaneous correlation function appearing in Eq. 
( 2.21| ), < u"[r(t)] u"[r(0)] cos 2 Oq r(t)r(0) > ro , with the initial displacement between the key solute 
and solvent atoms held fixed at tq = |ro| = 0.95a uv and 9$ = 0. As in Figs. |[-j3|, the reduced 
temperature ksT/e vv = 2.5. The solid and dotted lines both make use of the manifestly nonlinear 
solute-solvent Lennard-Jones potential u(r) in evaluating u"(r), but they take the time evolution 



to be given by the anharmonic dynamics of Eq. ( |2.22j ) and the harmonic dynamics of Eq. (3.1) 
respectively. The long dashed curve results from the corresponding linear INM theory: that is, from 
employing both the harmonic dynamics and the linear coupling assumption u"[r(t)] ~ it"[r(0)]. 
For Xe, a time interval of 0.1 tlj corresponds to about 350 fs. 
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FIG. 5. The instantaneous-normal-mode density of states for liquid Xe at a density 
p = 3.0g/ml and temperature T = 280K (computed by taking Xe to be a Lennard- Jones liq- 
uid at reduced density and temperature pa 3 = 0.95 and ksT/e = 1.26). Shown by the arrow, for 
comparison, is the vibrational frequency of I2 in its ground vibronic state. 
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FIG. 6. Predicted vibrational population relaxation times, T\, as a function of liquid density, 
p, for I2 dissolved in liquid Xe at 280 K. Results from the full instantaneous-pair-theory are 
indicated by dots with a line connecting them drawn to guide the eye. These values were obtained 
by averaging the instantaneous results over 100,000 liquid configurations for densities less than 2.0 
g/mL, 50,000 configurations for p = 2.2 g/mL, and 20,000 configurations for all higher densities; 
2 standard-deviation error bars are shown. The point indicated by squares are from the full 
molecular dynamics simulation results of Brown, Harris, and Tully (Ref. |7^) as reported by Harris, 
Smith, and Russell (Ref. |l0[). 
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